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We study a gravity theory where a scalar field with potential, beyond its minimal coupling, is also 
coupled through a non-minimal derivative coupling with the torsion scalar which is the teleparallcl 
equivalent of Einstein gravity. This theory provides second order equations of motion and we find 
large-distance non-perturbative static spherically symmetric four-dimensional solutions. Among 
them a general class of black hole solutions is found for some range of the parameters/integration 
constants with asymptotics of the form of hyperscaling violating Lifshitz spacetime with spherical 
horizon topology. Although the scalar field diverges at the horizon, its energy density and pressures 
are finite there. From the astrophysical point of view, this solution provides extra deflection of light 
compared to the Newtonian deflection. 


I. INTRODUCTION 

Teleparallelism M is an equivalent approach of General Relativity with deep implications which cannot be pro¬ 
vided by the standard metric formulation of gravity. For example, it provides a correct definition for the local 
energy-momentum tensor of the gravitational field which cannot be expressed in terms of the metric. This formula¬ 
tion is based on the use of the vierbein, along with a connection of vanishing curvature, as basic objects, not aiming to 
define another gravitational theory, but to describe anew Einstein’s gravity itself. Torsion turns out to be the appro¬ 
priate quantity for this reformulation. The Einstein-Hilbert Lagrangian is expressed (up to boundary quantities) in 
terms of torsion, so that only up to first order derivatives in the dynamical quantities (vierbein/connection) appear in 
the new Lagrangian. This Lagrangian is a particular quadratic combination of torsion which is called torsion scalar T. 
As a result, the Einstein equations have also their torsion analogue. As the metric formulation, also the teleparallel 
one is both diffeomorphism and local Lorentz invariant. By employing the Weitzenbock connection as solution of the 
connection equation of motion (vanishing curvature, but non-vanishing antisymmetric piece of the connection when 
expressed in coordinates, and hence non-vanishing torsion), the vierbein remains as the single dynamical variable and 
only the Lagrangian lacks local Lorentz invariance. Notice that in practice the teleparallelism condition of vanishing 
curvature is implemented through the Weitzenbock connection, and therefore the teleparallcl approach is different 
than other approaches (e.g. i) where torsion introduces a new dynamical field beyond the fields (metric, vierbein) 
describing the Einstein sector. Not only Einstein gravity but also Gauss-Bonnet gravity has been shown to possess 
its teleparallcl representation through another torsion scalar Tq 10]. Modifications of Einstein or Gauss-Bonnet 
gravity, different than the curvature based modifications f(R) [ll], [l2j], f(R,G) fl3l-fl5j have been constructed based 
on T, Tq, i.e. f(T) [l6l - [20j . f(T,Ta) [HJ theories have been investigated which in the covariant formulation are also 
diffeomorphism and local Lorentz invariant. 

Beyond the purely gravitational sector, there are in nature various sorts of matter fields. Scalar fields for example 
are predicted by fundamental theories and may be present in the universe at large scales mixed up with baryonic 
matter, and so be present in ordinary stars. Spherically symmetric configurations for self-gravitating scalar fields 
with various types of interactions have been derived as solutions of the coupled Einstein-scalar system S^-H^]. At 
cosmological scales the scalar field dynamics allows to investigate the features of the early universe, or a quintessence 
scalar field can be interpreted as the present dynamical dark energy sector j25j ] . The appearance of a non-minimal 
coupling of the form f(<p)R is motivated by many reasons, such as the variability of the fundamental constants, the 
Kaluza-Klein compactification scheme, the low-energy limit of superstring theory, e.tc. A non-minimal derivative 
coupling between matter and gravity is the next step of generalization, particularly if the Newton’s constant depends 
on the gravitational field source mass [26| . 

In [27| , a first step was made in an attempt to probe the effect of torsion to General Relativity when this is treated 
in its teleparallel reformulation, so the torsion scalar T was coupled to a scalar field <j> in four spacetime dimensions. 
Since for a diagonal vierbein, a simpler coupling of the form cj)T does not possess spherically symmetric solutions 
(however for three dimensions see |28|), a non-minimal derivative coupling of the form Tg^d^cfd^cj) was considered. 
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This is a novel coupling which provides only second order equations of motion. The metric-based analogous coupling 
of the form Rg^d^j) d„<j> gives higher order equations of motion [2(|, while the sector d v <j> of the general 

Horndeski Lagrangian also is a healthy theory with second order field equations [29;]. Of course, the above derivative 
coupling between torsion and matter is not the only one that can be constructed, but other combinations of d^d^cj) 


with quadratic (of even higher) powers of torsion could be presented. In 27] static spherically symmetric solutions 


were studied for the above coupling with torsion, however after the formulation was developed, only a general class of 
large-distance linearized solutions around its asymptotic Anti-de Sitter (AdS) form was found. Here, we will extend 
the analysis and find other (general and special) large-distance solutions in the non-linearized regime however, which 
turn out to be black hole solutions of the theory. The approximated solutions found have the meaning that they 
become of higher and higher accuracy as we restrict in particular regions of the parameter/integration constants 
space, which is not however fine-tuned or particularly narrow. The regime of applicability of these solutions is for any 
distance larger than some length scale determined by these parameters. A general such black hole solution turns out 
to have asymptotics which deviate from the AdS form and are of the hypescaling violating Lifshitz form. 

Our work is organized as follows. In Section II the basic elements and equations of the model appeared in [27] 
are briefly presented so that to make the present work autonomous. In Section III we proceed with the essential 
steps in order to simplify the master equation of the problem and facilitate our quest for finding large-distance non- 
perturbative spherically symmetric solutions of the theory. In Section IV we consider a special solution (i.e. for a 
particular value of the integration constant of the scalar field) and find black holes and other solutions along with the 
scalar field profile and the potential. In Section V a dynamical systems analysis reveals the global phase portrait of 
all the local solutions and serves to obtain an intuition for the approximations made in order to obtain the solutions 
of the next Section. In Section VI a general from the viewpoint of the number of integration constants black hole 
solution is found, its asymptotic behavior is analyzed and some astrophysical perspective of the solution is elaborated. 
Finally, Section VII is devoted to our conclusions. 


II. GENERAL FORMULATION 


We will consider the following non-minimal derivative coupling of a scalar field c/j with torsion 27 1 

1 


S = — —— / d 4 x e T — / d^x e 


G-eO 


g^d^d^ + v 


( 2 . 1 ) 


where V is the potential of the scalar field and the parameter £ has dimensions of length square, so ^/j^j introduces 
a new length scale. The torsion scalar T is defined by 


rj~< _ _ 


HvX 


j ^ _ rj~\ VflrjnX _ £1^LlsXrj n 


fii'X : 


where the tensor S^ x is 


g/ii/A _ LfcvXn _|_ p v _ gP v T pX ) = —S pX ‘ 


and 


K-flV A r) ^V[lX 


( 2 . 2 ) 


(2.3) 


(2.4) 


is the contortion tensor. The connection u> a b = uj a b ^dx p = w a bc e c is assumed to satisfy the teleparallelism condition 
of vanishing curvature R a bfJjU = ui a bll — ui a bfl v + u a cll uj c bll — co a cl ^oj c bfi = 0 (indices a, b, ... refer to tangent space, while 
/i, v,... are coordinate ones). In terms of an arbitrary orthonormal vierbein e a = (with dual e a = e a ^dx 11 and 

e = det(e a Al )), i.e. g^ = r] a b e a M where g a b = diag(—1,1,1,1) is the Minkowski metric, the torsion is defined by 
T a bc = uj a cb —oj a bc — e b /i e c ' y (e a Al jy —e a „ M ) and is a tensor under local Lorentz transformations and under diffeomorphisms. 
The metric is assumed to be compatible with the connection uj a bc , i.e. rf ab \ c = 0, and so u] a bc = ~bj baC i where | denotes 
covariant differentiation with respect to oj a bc . The importance of the torsion scalar T is that it provides after variation 
with respect to e a p the Einstein equations, i.e. eT equals eR up to a total divergence, where R is the Ricci scalar of 
the Christoffel connection T a ‘ i/A . 

A curvature-based analogue of the above non-minimal coupling is the term Rg^d^ipd^cf). The basic difference 
between these two is that while the curvature coupling provides higher than second derivatives (and therefore ghosts), 
here, the torsion coupling gives only second order equations of motion. In some sense, the coupling in (12.11) can be 
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said that it is closer to the curvature coupling G^d^d^cj) which also provides second order equations of motion, and 

indeed at the level only of background cosmology these two couplings coincide._ 

The equations of motion for the vierbein and the scalar field were found in [2?! to be 


(^-4^,,^) [(eS K At V),„eV + 


1 ^ 


4£ 


L 2 


eT^^+eSA^,^’*),, 


+ e^-cj ) ^»5 x - + VS X ) - ^-4^ tP ^yS dca uj bdc e a x e\ = 0 


(2.5) 


[e(l-2^T)^] /i - e — = 0, (2.6) 

where </> ,Al = g^ v <j) tl/ . As for the connection oj a bc , we adopt the Weitzenbock solution of the equations of motion 
Rabcd(u a bc ) = 0 , i.e. oj a bc = 0 in the frame e a M we will work with, therefore, w A = e a A e° v in all coordinate frames. 

We are interested in extracting static spherically symmetric solutions of the four-dimensional scalar-torsion gravity 
described by the action HU. Due to the spherical symmetry, the general form of the metric is diagonal 

ds 2 = - N(r) 2 dt 2 + K(r)~ 2 dr 2 + R(r) 2 dn 2 , (2.7) 

where dtt 2 = dd 2 + sin 2 ddip 2 is the two-dimensional sphere, and N{r ), K[r) and R(r) are three unknown functions. 
One could have chosen the gauge with R(r) = r, but it turns out to be more convenient to let the function R(r) free. 
As a next step in order to extract a vierbein that gives rise to the above metric, we choose the simplest diagonal 
vierbein of the form 


Substitution of 


1 

K 2 


e“ M = diag (N(r), I\{r) 1 ,R(r), R(r) sin#) . 

in the field equations (12.511 . (12.611 provides the following set of equations of motion 

2V\ _/ 1 ' 'W 
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( 2 . 8 ) 
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/(((/-C-") =0 ( 2 . 12 ) 


{ ,™v[i + ^-(- + aA) ]r -™.0, ( 2,3, 


NR 2 dV 


K d(f) 


where a prime denotes differentiation with respect to r. Equations (I2.9I) - (I2.12I) are respectively the non-vanishing 
(A, /z) = tt,rr, 90, dr components of the system (12.51) . An interesting characteristic, usual in teleparallel gravity, 
is the appearance of the off-diagonal equation (12.121) . although the metric ansatz (12.71) and the vielbein (12.81) are 
diagonal. This does not happen in the equations of motion derived from curvature-based actions, and in general, 
the extra equation imposes difficulties in finding solutions in the torsion formulation. Equations (I2.9I) - (I2.13I) have the 
reparametrizion invariance, so for r —> f(r) and K — > K^, N —> N, R —> R, </> —> <j>, the equations remain form 
invariant. This implies that one of these equations is not necessary because it arises from the others. Thus, we remain 
with four independent equations for four unknowns N(r), K(r ), </>(r), V (r) (i?(r) is not considered in the enumeration 
of the unknowns due to the choice of the radial gauge), so there is no arbitrary function left. As a result the previous 
system has basically a unique solution and this gives us a hope to extract some information on the structure and 
behaviour of the derived solutions not constrained by a prefixed potential. 

In [27j it was shown that if 


x = In R 


Y = 



(2.14) 

(2.15) 


where a dot denotes differentiation with respect to <f>, the system of equations (12.91) - (12.131) reduces to the master 
equation for Y ( x) 


2 — 

dx 2 


1 (dY\ 2 
Y V diJ 


+ 2 2 - 


77 v‘ 
~Y 


\dY o ■, 

■) — + 3 Y + 2gv 2 - 12+tT 
/ dx 


dY 

dx 


dY 

dx 


2 r ~2x 

— = 0, 


+ 3 Y + 2rjv 2 


(2.16) 
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where 


V 2^ 2 (1 2 (ft 2 !/ 2 ) 


^ i/ 2 (l — 6£ft 2 zA) 


(2.17) 


The integration constant v with dimensions of inverse length square is introduced through the relation <// = arising 
from (12.121) . After equation (12.161) is solved, the potential, the scalar field profile and the metric can be found. Namely, 
the potential V(R) is found from 


V = 


o~2x_^_±f 


1 /dY 


2rjv 2 


2 ij V dx 


3Y 


and the scalar field (f>(R) is obtained from 



The lapse function N(R) is found from 


where 


' dln(RN 2 )i 2 
dx 


z (x) 

Y(x) ’ 


Z = 


ArfY V dx 


j- + 3Y + 2rjv 2 


Finally the metric takes the form 


ds 2 = -N 2 dt 2 + 


dR 2 


/2R2 Y 


R 2 dttl. 


(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


Note that as far as £ and v remain unrelated, the coupling £ appears only through the combination £ari' 2 and this is 
due to the initial interaction term Tg^d^dy/j). 

It is also interesting to investigate the behavior of another physically meaningful quantity which is the energy- 
momentum tensor T“, of the scalar field <f>. This arises by variation with respect to e a M of the scalar field action, 

as 


i.e. of the second integral appearing in GED- One can find from (I2.9I) - (I2.12D the on-shell components of 
V t = T\ = \ ’ T R R =T r r = \ —^ ’ which express the energy density and the pressures of 

the scalar field. Since K(j)' = v, it arises T ( t = T e g = 2 > 


~t. t'Q 1 2V+v 2 J-R 'yr 1 2.V—v (l+8gi? ) 


6 £k ^-1 


III. SIMPLIFICATIONS 


We start by writing equation (I2.16|) as 


+ 2 ( 2 -^) — + 3 ( 1 -^ 

dx 2 Y V dx) V Y / dx V f? 2 


gv 2 \ dY 


Arj~ 


r niP‘ A- _— p 2x 

Y + 2 rju 2 + 24^ -- = 0 , 

y 2 §- + 3Y + 2 w 2 


which also takes the form 

;d 2 VY 


aVy- 


dx 2 


(2 

)- + 3fl- 

4ry 2 \ 

V Y j 

' dx \ 

fj 2 ) 


yv 2 + jA e ~ 2x 


dY 

dx 


+ 3F + 2gv 2 


= 0 . 


If we define 


u = Vy, 


equation m takes the form 

d 2 U / 0 yv 2 \dU 3/ _ V\ r/ V 3 y 2 yv 2 + ^e" 2 * 

dx 2 V U 2 ) dx 4 V fj 2 ) 2U fj 2 jLL Ijj _|_ 


(3.1) 


(3.2) 


(3.3) 


(3.4) 
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which is simpler than the initial equation (12.161) . 
Equation (13.41) for 


R » 


l VW\ 


becomes autonomous 


dx 2 


(2-< N 


K) 

V U 2 / 

1 dx 4V 

fj 2 ) 


V fj 2 ) 2U fj 2 dU q. Ijj q. lEl 


= 0. 


If 


dU 3 qv 2 

n = ^ + 2 U + ir 


(3.5) 


(3.6) 


(3.7) 


nfn 3 77 7 )v 2 \dn 1 2 3q 2 „ 3if 2 

n V } ~2 U ~~u)dU + 2^ + W V ~°- 


(3.8) 


equation (13.61) becomes 

' 2 ^° 1 ^ + ^ 

7^ 77^ 

It is convenient to define a = qv 2 , /3 = > 0. The parameters a, /3 are not independent, but they obey the relations 

a = 2 (i- 2 2 ") ’ P = l( i- 2 ^ ) 2 = Kl^ -3 ) 2 ’ w here 5 = ^n 2 v 2 . It is seen that for a > 0, the parameter /3 can take any 
positive value, while for a < 0 it is /3 > The inequality (13.51) gets the form 




'l\/M 


= \ 2 


H- 


Finally, equation (13.81) becomes 


Equation (13.101) is also written as 


or also 


dtt 1 


i >{ n -2 U -ij)du + 2 a ‘- l3un + al3 = 0 ' 




dn _ 

du ~ n -1 u-§ 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


Although (13.121) is a first order differential equation, significantly and unexpectedly simpler than the initial equation 
(12.161) . there is no known method how to solve it. However, its characteristic form with the presence of the inverse 
powers of U, will enable us to find “large”-distance non-perturbative solutions. As always, the word large is meant in 
comparison with scales and integration constants of the problem, and depending on them, the corresponding distances 
can be relevant for a physical situation. 

Note that the system of equations (I2.9I) - (I2.13I) for vanishing scalar field <j> = 0 and a cosmological constant as the 
potential V has the standard Schwarzschild-(A)dS solution. A vanishing scalar field means that v = 0 and then it is 
a = P = f, | while Y, U of (13.21) . (13.41) become infinite. Moreover, the autonomous equation (13.61) cannot 
recover this limit since the inequality (13.91) is never satisfied. The reason is that the exponential factor in (13.41) becomes 
significant in this limit. However, although the general solutions to be derived will not be continuous deformation of 
Schwarzschild, this does not necessarily mean that they will not capture some interesting phenomenology, for example 
they could be related to galactic scales or beyond. 
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IV. SPECIAL SOLUTIONS WITH v ->■ 0 


In the limit v —> 0 the theory posseses solutions, including black holes, with non trivial scalar field and potential. 
For v —> 0, the scalar field <j> 0 and (12.191) is not defined due to that dcj> vanishes. However, the “normalized” field 

<j> = ^ obeys 4>' = K which is a meaningful equation. Similarly, although U —> oo in (Em), however this equation 
still makes sense. Indeed, if we define the “normalized” variable U = vU, (13.41) becomes in this limit 




V dx 


-U 


d 2 U 


dx 2 


+ 2 


<UA' 


dx 1 


+ 4 6 


= 0 , 


where there is no free parameter left in (14.11) . Then, (12.191) becomes 


(4.1) 


i = ±u 

dtp 


and equations (12.201) . (12.211) are meaningful. Finally, the metric (|2.22D becomes 


(4.2) 


ds 2 = -N 2 dt 2 + 'F n + R 2 dfll. 
FPU 2 

The potential V is found from (12.181) to be 

v =U e ~"-™ 2 - 2 u t)- 

We will find as a warm up large-distance solutions in the above limit. 

For F = e x U we get from S3 


If v(F) 


dF 
dx 5 


then 


or equivalently 




V dx 


dx 2 


F 




/ l w 

/ dv 

1 

( v 4— F) 

(v—— — F ) 

H—~ 

V 2 / 

V dF ) 

4 


dv & 1 

v —— = r -. 

dF 2jP+4u 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


The variables F, v are dimensionless and we can check by numerically solving (TO) that every solution extends to 
large values of |F|, |u|. Thus, for |F|, |u| >lwe get from (14.71) 


with general solution 



= F 2 — a , 


(4.8) 


(4.9) 


where a is integration constant. Then, 


F = 


c tR 2 +t 2 
2tR 


or 


F = 


r 2 R: 2 + a 
2 tR 


(4.10) 


where r is another integration constant. Depending on a, r, this solution can indeed correspond to large values of R, 
since we see that large values of R can correspond to large values of |F|, |u|. Integrating equation (12.201) we find N. 
Finally, rescaling the time t, there arise two sorts of metrics 


ds 2 = — F 2 dt 2 + + R 2 dFl 

F 1 


2 

2 > 


(4.11) 
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or 


ds 2 


1 

R 2 F 2 


dt 2 + 


dR 2 

~pT 


R 2 dfli 


(4.12) 


Thus, there exist four different metrics. 

Depending on the constants a, r, the metric (14.111) may possess horizon. To be concrete, let us consider the first case 
of (14.101) with r > 0, a < 0 and \<r\ -C 1. Then, it is obvious that large values of R correspond to F -C 1 (or |F| 1), 

which means R From (14.91) it is also |u| 1, while the quantity |2F+4u| is large even for Fv < 0. On the other- 

hand, the horizon of (14.111) corresponds to F = 0, which means i?hor = J ■ This value can be made arbitrarily 

V M 

larger than ^ by choosing \a\ sufficiently small, and so the existence of the horizon within our approximation has 
been shown. In this case, (14.111) with the first of (14.101) is the black hole 


„ (|cr|i? 2 — t 2 ) „ 

ds 2 = - Vl ' ; dt 2 


At 2 R 2 


T dR 2 + R 2 dnl 


^ 2 R 2 (| \a\R 2 -r 2 Y 

and the solution is certainly valid outside the horizon. Asymptotically, it gets an AdS form 

,9 er 2 9 4r 2 dR 2 

ds °° = dt " + 72‘ _ R2 _ + R dfl 2- 


(4.13) 


(4.14) 


The length scale introduced here is £ e g = -fpz and the corresponding effective cosmological constant A e ff = —3£ e g. 
Similarly, for the second case of (14.101) with r < 0,<r < 0 and \a\ -C 1, it is obvious that large values of R correspond 
to F -C 1 (or \F\ 1), which means R S> From (14.91) it is also |u| 1, while the quantity |2F+4u| is large even 

for Fv < 0. On the other hand, the horizon of (14.111) corresponds to F = 0, which means i?hor = -pjrp. This value 
again can be made arbitrarily larger than for |cr| sufficiently small, and so in this case (14.111) with the second of 
(14.101) forms another black hole 


ds 2 


(r 2 R 2 — |cr|) 2 
4r 2 i? 2 


dt 2 + 


4 t 2 R 2 

(r 2 i ? 2 —|cr |) 2 


dR 2 + R 2 dn% . 


Asymptotically, it also gets another AdS form 


(4.15) 


ds 


2 

oo 



R 2 dt 2 + 


4 dR 2 

T 2 ~EP 


+ R 2 dn% 


(4.16) 


with ej = £. 

Concerning the scalar field <fi, we find from (14.21) for the two black hole solutions (14.131) . (|4.15l) respectively 

^ — <^o = ±—In (|cr|7? 2 —r 2 ) or (j — cj)Q = ±— In (t 2 R 2 — |cr|). (4-17) 


Notice that the scalar field diverges at the position of the horizon and also asymptotically. The scalar field is a 
secondary hair since it does not introduce a new non-trivial integration constant (i.e. <f> is not varied independently 
of the black hole parameters). 

Finally from (14.41) . since we find the potential V respectively 


g 2 F~ 2 —l+2e^/l + |g|F 1 ~ 2 or t 2 F- 2 -l+2e v /l + |g|F 1 ~ 2 

k2t ' 2 (1 + Vl + kl F~ 2 ) 2 OT k2 (l + y/l + \a\F~ 2 ) 2 


(4.18) 


where e is a ± sign. From (14.171) . converting R to <f>, we can easily find F((f>) respectively 

4r 2 r 


F~ 2 = 


e =Ff +T 2 e T^-(<l>-<l>o) 


or F ~ 2 = At 2 


g=FT(0-0o) _ cre =F2T(0-0o) 


(4.19) 


Substituting (14.191) into (14.181) we find V(</)). Note that for any solution with specific a, r the corresponding potential 
V that supports the solution also depends on a, r. We could rescale to (j) = ^ (<p — </>o) or </) = r(0 — <j) 0 ), but still 


















































FIG. 1: Phase portrait (U, Cl) for a = — |, j3 = 27 


a, t remain in V. At spatial infinity it is F 2 —> 0, therefore the potential becomes for the two cases V —> 


(2e—l)fT 2 

4 k 2 t 2 


or 


V — > ^ 2e ^ 2 T , and for e = — 1 these values coincide with the corresponding values of A e g/n 2 . For e = 1, these values 
of V are positive, therefore the cosmological constant coming from the potential in the action is positive, while the 
effective cosmological constant A e g of asymptotically AdS space is negative. This discrepancy is due to the growth of a 
non trivial profile of the scalar field which is coupled to the torsion and modifies the asymptotic form of the spacetime. 
In [27} a similar discrepancy between A e ff and V has been found, but not for positive V which is the case here. A 
positive value of V has the merit that it can be related to the vacuum energy of spacetime. On the other hand, in 
1231 and other solutions the asymptotic value of V is the same as A e g/ k 2 . Note also that at the horizon the potential 

V is finite. Concerning the components of the energy-momentum tensor we have T t t = = l~ T r = T e g = —V. 

Therefore, although the scalar field diverges at the horizon, its energy density and the pressures are finite there. 
Asymptotically, still the energy density of the scalar field and the pressures remain finite. Finally, notice that due to 
the higher order pole at the horizon, the temperature vanishes in analogy to the Reissner-Nordstrom solution. 


V. DYNAMICAL SYSTEMS ANALYSIS 


The second order autonomous differential equation (13.61) can be converted into a two-dimensional dynamical system 
with one function the variable U and as second function some combination containing the derivative e.g f l. To 
study the system it is better to convert into dimensionless variables. We define 



K 



K 



(5.1) 


where U > 0 and Cl is real. So, instead of the equation (13.61) we have the following equivalent system 


dU 

dx 


= Ci--u- “ 

2 U 


dCl 

dx 


pu'n-W 
2 n 


(5.2) 

(5.3) 


where /3 = |(4d — 3) 2 . We can find the fixed points of this system by setting ^7 = = 0. Then, we find these fixed 

points to be 

= /±2^3?^ ’ n * = ± \j ±2VP- 3 ■ (5 ' 4) 

For a < 0, /3 > ^ there is exactly one fixed point (17* = ^ 2 ^ 3 + 3 ’ = “ \J There are two negative 

eigenvalues 3Ai_2 = — 2(^/3/? + 3) ± \J 12/3 + 6-v/3/3 + 9 of the linearized system, thus the fixed point is attractor 
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FIG. 2: Phase portrait (U, Cl) for a =a /3 = 21 



FIG. 3: Phase portrait (U, Cl) for a = |, j3 = A, 

(stable node), as also seen in Fig. 1. In Fig. 1, 2, 3 the horizontal axis is U and the vertical Cl, and the arrows 
show increasement of the radius R. Moreover, due to the pole Cl = 0 in (15. 31) . the phase portraits are separated 
into two quadrants, the upper one with Cl > 0 and the lower with D < 0. For a £ (0, i) U (l,+oo), /? > | 

there is again exactly one fixed point (17, = yj 2 ^|g„ 3 ", = \j case > the two eigenvalues 3 Ai j2 = 

2(\/3/3 — 3) ± \J 12/3 — 6\/3/3 + 9 have opposite signs, therefore the fixed point is saddle, as also seen in Fig. 2. Finally, 
for_4_< a < 1, 0</3<| there is no fixed point, as it is seen in Fig. 3. The values (15.41) of U are the same found 
in (27| studying the asymptotic behaviour of equation (12.161) . From the above analysis it becomes clear that only for 
a < 0 (case A of [27j | j t he fixed point is attractor and indeed corresponds to a large-distance asymptotic solution. 
For a > 0 (case B of [23]) the fixed point is saddle and does not correspond to an asymptotic solution. It would be 
interesting for a < 0 to investigate beyond the linearized asymptotic behavior found in [27| , also the non-perturbative 
regime of the solution, however, this will not be studied in the present paper. In the present paper we will find the 
large-distance non-perturbative solutions which correspond to large values of U, as shown in the upper quadrants 
of Fig. 1 and 2. These solutions will be now shown that are attracted by stable fixed points at infinity. 

Due to that the dynamical system (15.21) . (15.31) is non-compact, there could be non-trivial fixed points at infinity, i.e. 
at the asymptotic region of the phase portrait ( U,Cl ). These can be studied using the Poincare projection method. 
We define the new variables (coordinates of phase space) by 

U = ———^ cos 9 , Cl = -- sin 9 

1—r 1 —r 


(5.5) 
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with — j <&<j,0<r<l. The upper quadrants of Fig. 1, 2, 3 correspond to 9 > 0, while the lower ones to 9 < 0. 
The limit f —> 1 _ corresponds to infinite distance in phase space, U 2 + fi 2 —> oo. In terms of r, 6 the dynamical 
system (15721) . (ESI) becomes 


dr 

dx 

d6_ 

dx 


r— 1 

1 ■ 
2f2 . 


|2d/3(l —r) 2 + 2d(l —2f) + f 2 [2(l + o;)+cos20—(/3+1) sin2d] j 
2d(l —f) 2 (tan 0—/3 cot 0) + ((3+l)r 2 cos 2 8 + f 2 (sin 29 +/ 3 — 1 ) . 


(5.6) 

(5.7) 


As r ~ 1 the leading terms of equation (15.71) for are 


dd 

dx 


1 

2 


[(/3+1) cos2d+sin20+/3— l] , 


(5.8) 


while at linear order 4^ = 0. The critical points 0* at infinity are obtained by setting (4=0 in (15.81) and solving for 
0, thus 


(/3+1) cos20*+sin20*+/3—1 = 0 . (5.9) 

This equation has for any fd > 0 two roots for d*, one positive and one negative. Therefore, there are always two 
fixed points at infinity, one in the upper quadrant of the (r, 9) plane and the other in the lower quadrant, as are also 
depicted in the uncompactified plots of Fig. 1, 2, 3. Since = 0 we cannot rely on the linearized analysis to 

examine the stability of these fixed points and numerical examination is needed. Indeed, it can be seen numerically 
that for the parameters of Fig. 1, 2 a stable fixed point exists for 0* > 0 (for which it is -4((=) < 0) and an unstable 
fixed point exists for 0* < 0 (for which it is -4(44) > 0), in agreement with the phase portraits of Fig. 1, 2. 

du 


VI. GENERAL SOLUTIONS 


Equation (13.121) can be approximated in the case that U,tt > 0 are large (compared to the gravity coupling k), or 
more precisely the dimensionless quantities U, O > 0 are large. It is obvious from Fig. 1, 2 that large values of U, II 
correspond to large R, something that will be verified after the solution is found. More precisely, for 


we get 


^ 3 rr 

|a| 

1 ^ 

n - -u 

> rr ■ 


2 

U 

2 (3 


dfl _ 2/3U - O 
dU ~ 20 - 3U ’ 


( 6 . 1 ) 

( 6 . 2 ) 


which is a homogeneous equation. The above two inequalities are obviously consistent for large U , O, and moreover, 
these inequalities (together with (13.91) 1 will define the exact i?-domain where the solution is valid. The puzzling 
situations with O ~ |(7 or 9 ~ 2(3U (U large) do not occur since they provide for large R that ^ oo or 0 
respectively, which is seen from Fig. 1, 2 not to be the case. The above inequalities set restrictions in the allowed region 
of the 3-dimensional space (x, U, ^4) of the initial second order differential equation (13.61) where the orbits reside. Due 
to this, although we will find a general large-distance (non-linearized) solution with the correct (maximum) number 
of integration constants, the solution will fail to describe other regions in the space of initial data which also provide 
large-distance solutions. This is obvious from the dynamical systems analysis of the previous section. 

Assuming, for example, the parameter \a\ to be sufficiently small, which means |S| = |£|« 2 z/ 2 ^ 1, the inequalities 
can more easily be satisfied and the validity of the approximation becomes even more extended. At the same time 
the value of the parameter /3 is very close to 27/4. Such values |S| 1 for the parameters were argued in [27] that 
could in principle reduce a large value of the vacuum energy to a small effective cosmological constant. Moreover, for 
such S the inequality (13.91) becomes R 2^/[|[, so this inequality is determined by the length scale defined by the 
non-minimal coupling £. Additionally, since v is related to the integration constant of the scalar field <j >, it is expected 
to take large values for macroscopic solutions, thus even for small |£| it can be |E!| ^ 1. However, we do not restrict 
our analysis only in this range of parameters. 
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Setting 


equation (16.21) becomes 


which is separable. Setting 


(with $ > — 5 ), equation (16.41) becomes 


^ n „ 
$=- >o . 


d<k 

u w +2 


<h 2 -$-/3 
2<k —3 


= 0 , 



d’k 4» 2 -7 2 
U dU + "1>-1 


= 0 , 


(6.3) 


(6.4) 


(6.5) 


( 6 . 6 ) 


where 7 = P+\ > 5 (and thus ^k > 

R) takes the form 


— 7 ). Equation (13.71) which makes the connection with the radial variable x (or 


<i'k 4> 2 —7 2 / a \ 

= \k —1 ( 1_ + U 2 ) ' 


(6.7) 


The parameter a, which had temporarily disappeared in (16.2D . came up here again. The quantity U 2 in (16.71) will be 
found as a function of W by integrating (16.61) . Indeed, the solution of (16.61) is 

U 2 = C (T+y ) - ” -1 |\k — 7 |~ -1 , ( 6 - 8 ) 


where C > 0 is integration constant that distinguishes the solutions in the (C7, fl) space. In terms of f l the first integral 
(16.81) takes the form 



(6.9) 


From (16.81) it is seen that only for 'k ~ 7 > 1 it is U —> 00 , in agreement with Fig. 1, 2. Therefore, we focus on 
7 > 1 (which means @ > |, while a < |- or a > k 2 ) where U can indeed extend to infinity. Furthermore, since 
II = (4/ + \)U, for such 7 there are solutions where both U, f l can extend to infinity. From (16.81) it is seen that there 
are two classes of solutions, one with 'k > 7 and the other with < 7 (these solutions are represented respectively 
by the left/right orbits of the upper quadrants in Fig. 1, 2). From (16.61) it is seen that for the solutions with 4* > 7 
it is T a monotonically decreasing function of U. To see if the above solutions correspond to large i?, we set \k ~ 7 
in the integral (16.81) or even the equation (16.71) . Indeed, we get R ~ i?o/|'k — 7 |*r, so the region R —> 00 in included 
and the metric can be extended to large distances. Obviously, <k is not necessarily close to the asymptotic value 7 in 
the range of applicability of the solution and the constraints (13.51) . (16.11) will determine the exact domain of R. 

The metric (12.221) becomes 


ds 2 = -N 2 dt 2 + 


|\k — 7 | 1 7 ('k+ 7) 1+7 dR 2 


2 CR 2 


R 2 dVL 2 2 , 


where the relation of R, \k in (16.101) is found from (16.71) 


i?.(\k) = R 0 1 -*+»r 2 = R 0 e " 


( 6 . 10 ) 


( 6 . 11 ) 


with Rq > 0 integration constant. Integration (16.111) cannot be performed analytically for arbitrary 7 . However, we 
will see that this is not necessary because within our approximation, equation (16.111) becomes simplified. Metric (16.101) 
has an implicit form through integrating and inverting (16.111) for *k(i?). However, we can obtain the metric (16.101) in 
a more explicit form without the need for inversion (although not in the standard radial gauge) as follows 


ds 2 = -N 2 dt 2 + 


('k-lfd’k 2 

v 2 C |\k — 7| 1+7 ('k+ 7 ) 1- ~ (1 — 'k + a/7 -2 ) 2 


i?(>k) 2 ^2 ■ 


( 6 . 12 ) 
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Since the master equation (12.1611 is of second order, there are two extra integration constants C, Rq in the solution 
(16.111) . beyond the integration constant v. Therefore, after the lapse function N is found, the metric (16.101) or (16.121) 
will be the general approximate solution in the domain of its validity. 

The lapse metric function N can be obtained from the equations (12.211) . (12.201) from where we obtain 


d In (RN 2 ) efj 1 fdY \ 

—5^ = 2?y(^ +3y+2 “) 


with e = ±1 a sign symbol. Integration of (16.131) gives 




N = - R 3 U 


>3rr2 2a 


where c > 0 is integration constant. Using (ED , liuTD - m we find 


N 2 {R) = 


R 3 e J 


_ £ 2 . 
2 V 


R Ljvlr-^l 1 7 7 - 


(6.13) 


(6.14) 


(6.15) 


where 


J = 2a 


\F-1 d\ 1/ „ f \F-1 d\ 1/ 

v [/ 2_ 7 2 a +( i _^) j72 " / T 2- 7 2 j + £( i _ W )| W - 7 | 4- 1 ($-|- 7 )-^- 1 


(6.16) 


~ £ 2 . 

and c = cC 2 ’ >0isa redefined integration constant. The integration constant c is not essential since it can be 
absorbed into a redefinition of the time coordinate t. 

The potential can be found from (12.181) . after use of (13.71) . (16.31) . to be 


i ' v, =-(® + \y + h J. 

C(T + i) 


T- 7 I 1 r ('F+y) 


2 v 2 R 2 
a 
2 


1 


l + V 


2 v 2 R 2 ' 


Finally, the scalar field configuration can be determined from (12.191) using again (13.71) . (16.31) 


= </>i + ei 




tf -1 


d>F 


T 2 - 7 2 l-'F + at /" 2 U 

r T —1 |i£_ 7 |!( 1 -y) (^ + 7)5(!+y) 


\ff2_~2 


T l-4' + §|'F- 7 | 1_ “ ($+ 7 ) 


Y+i 


■ d <F, 


where is another ± sign and <pi is an integration constant. 
In terms of 4/ the inequalities ED are written as 


(T+7) 1+7 1 |T —7| 1-7 " 

I'F-H 


„ c_ ( 4 /+ 7 ) 1+7 ~ IT- 7 I 1 - 7 - C_ 

H ’ P (T + i)|T-(27 2 -l)| \a\ 


(6.17) 

(6.18) 


(6.19) 

( 6 . 20 ) 


( 6 . 21 ) 


From the second condition of (16.211) it is seen that the branch with *F > 7 has a pole at *F = 2y 2 — 1 > 7 . Therefore 
choosing a value for ^ both conditions (16.211) will be satisfied for a solution which is defined from \F = 7 (for R —> 00 ) 
up to some value smaller than 27 s —1 (that corresponds to the minimum R). If ^ is sufficiently large, \F can approach 
the value 2y 2 — 1. Similarly, from the first condition of (16.211) the branch with \F < 7 has a pole at 4/ = 1. Therefore 
choosing a value for ^ both conditions (16.211) will be satisfied for a solution which is defined from some value larger 
than 1 (that corresponds to the minimum R) up to \F = 7 (for R —> 00 ), since from (16.61) *F is now an increasing 
function of U. If ^ is sufficiently large, \F can approach the value 1. 

Hopefully, it is obvious that the first of the conditions in (16.211) assures that the quantity ^(’F+y ) 1+7 1 1— y | 1 1 
in the integral (16.111) is negligible. Then, the integration can be done and it gives 


R — Rq 


*F+7 


*F — 7 


( 6 . 22 ) 
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Accordingly from equation (16.81) 


2 7 \RqJ 




(6.23) 


where s = sgn(\I/ — 7 ) discerns the two branches with $ > 7 or $ < 7 . In the branch with d' > 7 , it is seen from 
(16.221) that as T increases, the radius R decreases, while in the branch with d' < 7 , as ’I' increases, R also increases, 
in accordance with the previous analysis. 

Finally, the conditions (16.211) take the following form in terms of R 



(6.24) 

(6.25) 


respectively. For s = 1, these inequalities say that R can become as close to Rq as we wish, given that /3|a |/C is 
sufficiently small. Since /3 > | this last condition means \a\/C <C 1 and the inequalities (16.11) are satisfied practically 
for any R larger than Rq. Furthermore, if Rq |i'| _ 1 |a|~ 1 / 2 , the condition (13.91) is satisfied for any R outside Rq. 
Therefore, all the conditions have been satisfied. For example, if |£| \Pv 2 1, the condition \a\/C <C 1 becomes 

4C|£|^ 2 1 which is satisfied if C is not particularly small; moreover, for Rq 2^/j£|, the condition (13.91) is satisfied. 

To proceed with the solution, now the integral J in (16.161) is found to be J oc ^ f ( |^ 7 | )^d l F, which around the 

dangerous for divergence point T = 7 gives J oc — 7 I 1-7 . Since 7 > 1, it is |\I/ — 7| 1-7 finite and it will be 

J ~ 0 for \a\/C 1. Using (16.221) . the metric (16.101) becomes 


ds 2 = — 


(~ 

\Rq 




1 — 5 


Ro\ 2 ^ 


/£to 

V R 


£\/f 


dt 2 + 


dR 2 


CR 2 


where ( = esgn[(l — 22)(1 — 62)] = ± 1 , C = yy/ > 0 and t = 


(ff 7 


t. 


R 2 dn 


2 > 


(6.26) 


-sC \[Jii 2 


Rn 


We will be interested in the solution with s = 1 (*F > 7 ) which presents a pole. As explained above this solution is 
valid for distances immediately outside Rq and for this reason, Rq for ( = 1 will be called horizon. The metric (16.261) 
in this case takes the form 


ds 2 = — 


(R 


1- 

fRo\ 2 ~<' 

Vi?o 

) 

V R ) 


dt 2 + 


dR 2 


CR 2 


(if [l-(f) 2 f 


R 2 ddl. 


(6.27) 


where for convenience we remind that C > 0, 7 = \J~P+\ > 1 and /3 = f(^^3) 2 >| with a < y or a > k 2 (the 

parameter a is related to the coupling £ and the scalar field integration constant v through the combination ^k 2 v 2 ). 
The analysis of the curvature invariants reveals that these are all finite at the horizon Rq and diverge at infinity. The 
integration constant Rq is expected to be related with the mass of the black hole. Asymptotically the metric (16.271) 
takes the form 


ds 2 = - 


Vi?o 


V 


3(2-y+1) 

27 — I 


dt 2 + 


CR 2 


dR 2 
' JL\ 21 

< Rq / 


+ R 2 dnZ. 


(6.28) 


This metric defines an asymptotic behaviour different than that of AdS space. Similarly to AdS, the lapse function also 
goes to infinity for large distances, but here the scaling behavior is different. More precisely, note that the exponent 

2y is larger than 2 and approaches the value 2 as 7 approaches 1. On the contrary, the exponent \J 3 — 1 
is smaller than 2 and positive, and approaches also the value 2 as 7 approaches 1. Therefore, for 7 very close to 
1 the asymptotic metric (16.281) gets close to AdS, while as 7 departs from the value 1 , the metric gets a different 
structure. For example, for 7 = 7/2, the lapse function in (16.281) becomes proportional to R (to give an astrophysical 
perspective, the metric (16.281) seems to provide extra attraction at large distances, while fittings of linear potentials for 
exponential galactic disks have been shown to explain the almost flat galactic rotation curves [30] and such potentials 
yield galactic stability without the need of dark matter (3j|). For a circular orbit at constant radius R the conditions 
EN 2 = 1 — jyjyy, j 2 = ^ 37 T have to be satisfied Q, where the constant E is the energy of the moving particle 
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(E > 0 for material particles, E = 0 for photons) and the constant j is related to the angular momentum. For the 
metric (16.281) both conditions are satisfied since 7 > 1 , therefore at large distances circular orbits are supported. 
Making the transformation 


p = R\ 


the metric (16.271) takes the form 


ds 2 = p d 


where 


i?n 7 \Vf ,.2 , R 0 


'/->) 


dr z 


dp 2 


+ -7T!-— £ W + P ^ 2 

Cl 2 p*( l -^) 2 


(l--) < 2 

0 < 2 = 1 - J - 

3-1 

/ 3(2 7 -|-l)" 

V 7 / 

2 7 


v 27-1 


< 1 


(6.29) 


(6.30) 


(6.31) 


and t = Rq * 27-1 t. The metric (16.301) is a hypescaling violating black hole [33| with spherical horizon topology. 
A hyperscaling violating black hole is a generalization of the Lifshitz black hole where d = 0 (in our solution the 
metric can never asymptote a Lifshitz spacetime). The Lifshitz metric arises as solution of gravity theories with 
negative cosmological constant coupled to appropriate matter with the simplest such theory also including an abelian 
gauge field [34) (a pure Einstein gravity with cosmological constant cannot produce an anisotropy in spacetime). Such 
metrics have also been found as solutions in string theory and supergravities which arise from string constructions 
[35| . Effective gravity theories with a Maxwell as well as a dilaton field (in general a scalar field with a non-trivial 
potential) are quite rich and have been shown to contain hyperscaling violating solutions [dil . Note that the majority 
of the Lifshitz or hyperscaling violating Lifshitz solutions and the corresponding black holes in the literature have 
planar (horizon) topology and are assumed to have direct correspondence with condensed matter physics through the 
AdS/CFT conjecture (for a spherical horizon topology to be obtained an extra gauge field should be added). On 
the contrary, the spherical symmetry found here may offer to the solution some significance at local astrophysical 
objects at large distance scales. In general it is not easy for a given theory to possess Lifshitz or hyperscaling violating 
solutions, and as referred, for example the introduction of some extra matter source or higher order gravity theories 
are required. Here our scalar-torsion theory is an additional case which provides such solutions and notably these 
solutions are general, while in most of the known cases the arising such black holes are special solutions. According 
to the standard notation, d is the hyperscaling violation exponent, while z is the dynamical critical exponent which 
indicates the anisotropy between time and space. For our solution the values that these parameters can take are seen 
from the conditions (16.311) . The asymptotic form of (16.301) is 


= P 


- 1 ? 


- p 2z df 2 + 


Rf^ dp 2 

cy p 2 


+ p 2 dfll 


(6.32) 


The scaling transformation t — > X z r,p —> X~ r p,Xi —> Xxi does not act as an isometry for the metric (16.321) . so (16.321) 
is not scale invariant, but it transforms conformally as ds ^ —> X^ds"^. 

In the context of AdS/CFT, a non-vanishing d indicates a hyperscaling violation in the dual field theory. In the 
four-dimensional framework we are working, theories with hyperscaling at finite temperature have an entropy density 
which scales with temperature as S ~ T*. For hyperscaling violation there is a modified relationship S ~ T ~^~, 
indicating that the system lives in an effective dimension d e e = 2 — d. For the present solution it is 0 < d e s = ^ < 2. 

The potential accompanying the black hole (16.271) is found from (16.181) to be 


V = 


/ R \2(7-l) 

1 - 

fRo) 2 ^ 




\r) 



' Ro\ 27 


v 

T 


1 


(6.33) 


2 aR 2 ' 

It is seen that at the horizon the potential is finite, while at infinity it diverges as V ~ (R/Ro) 2 ^~ 1 ' 1 . Finally, the 
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scalar field associated with the metric (16.271) is found 1 to be 


eiM / .Ro y- 1 

( 7—1 )Vc^ R ' 



(6.34) 


Since the parameter 7 depends on the integration constant v of the scalar field, the scalar field (16.341) is a primary 
hair. The scalar field diverges at the horizon and is finite at infinity. This behaviour at infinity is different than 
the behaviour found in [271] for AdS asymptotics, where the scalar field evolves logarithmically with distance. At 


infinite distance the scalar field behaves to dominant order as <j) — </>i ~ — 
corresponding behaviour of the potential V ( cf >) for 4> — (j)\ ~ 0 to be 


eiM 


(1-llV^ R 


Aayy-1 




Then, we can find the 


v{4>) * - 


27 +I v 2 
2 ( 7 -l ) 2 a((f>—(f>i) 2 ' 


(6.35) 


Thus, the potential V{<j>) close to the origin <j> — <j>i = 0 is very steep and can be either positive or negative depending 
on the value of a. At the opposite limit of distances close to the horizon, the potential gets an almost constant value 

Concerning the components of the energy-momentum tensor , although the scalar field diverges at the horizon, 
it is obvious that its energy density T f t and the pressures T~ R R , tg are finite there. Asymptotically these components 
of diverge. 

To summarize, the most significant solution found in this section is described by the metric forms J6'.£7[ ), \6.30\) 
with the asymptotic structures H6.28\) . \6.32\) . The solution develops an horizon Rq and R is practically defined outside 
Rq when Rq |^| -1 |a| -1 / 2 and -C 1. The solution is supported by the potential \6.33\) with the scalar field 


profile \6.3f\ ). 

The Hawking temperature T of the black hole is determined by the periodicity of the Euclidean metric ds 2 E = 
g TT dr 2 + g EE dR 2 + R 2 dH, 2 = N 2 dr 2 + I\~ 2 dR 2 + R 2 dfl\ obtained by the analytic continuation t = —it. Thus, T is 
given by the standard formulae 47 tT = ( dg TT /dR )/ y/g TT gRR\ r 0 = sjK 2 /TV 2 d(N 2 )/dR\n 0 = 1 / N 2 /K 2 d{K 2 )/dR\n 0 = 
yJ[d(N 2 )/dR\ [d(K 2 )/dR\\n 0 , given that K 2 vanishes at the horizon R 0 . Due to the higher order pole at the horizon, 
the temperature vanishes in analogy to the extremal Reissner-Nordstrom solution. 

From the astrophysical point of view, the motion of a freely falling photon in the static isotropic gravitational field 
(12.71) is described by the equation (dJd) 2 = R A K 2 ( — Ay). Since the field is isotropic, the orbit of the particle 

can be considered to be confined to the equatorial plane 9 = tL . At the distance R* of closest approach to the center 
it is ^ = 0, thus j 2 N(R *) 2 = R 2 . The deflection of the orbit from the direction of initial incidence at infinite 
distance is A ip = 2| ip(R*) — y>oo| — 7 r, where tpoo indicates the incident direction. The larger the quantity (^) 2 in 
the previous differential equation, the larger the deflection angle is. It can be easily seen that at large distances 
(^A ) 2 for the metric (16.271) is bigger than that of the Schwarzschild metric. Thus, there is for our solution an extra 
deflection of light compared to the Newtonian deflection. The situation of increased deflection compared to that 
caused by the luminous matter has been well observed in galaxies or clusters of galaxies. In general, the metric 
ds 2 = —R 2a dt 2 + g R 2 b + R 2 dTl\ = p~' } (—p 2z dt 2 + + p 2 dtd |), d = 2 — |, z = 1 + {p = R b ) has an extra 

deflection of light at large distances if a < 1,6 — a>0 which mean (2 — $)(1 — z) > 0, (2 — d)(2 — 2z + d) > 0. These 
conditions are obviously satisfied for the solution found here. 


VII. CONCLUSIONS 

In this work we have extended our previous analysis [27 ] on the quest of finding spherically symmetric solutions of 
a scalar-torsion theory. More precisely, we treat the torsion not as an independent field but the teleparallel condition 
is imposed as a constraint, therefore the corresponding connection is assumed to have vanishing curvature. This last 


From the Gauss’ recursion formula c(c + 1) 2 F 1 (a, b', c; z) — c(c + 1) 2^1 (a, b; c + 1; 2 ) — abz 2 -Fi(a + 1, b + 1; c + 2; z) = 0 (p. 1010 of 13711 ) 
and 2 F 1 (a, b; a; z) = (1 — z)~ b [3Sll . we get the relation (a + 1 ) 2 F 1 (a, 6 ; a + 1 ; z) + bz 2 F 1 (a + 1, b + 1 ; a + 2 ; z) = (a + 1)(1 — z)~ b . Since 


^ 2 Fi(a,b-,c;z) = ^ 2 F 1 (a + 1, 6 + 1; c + 1; z) 


we get 4z\z a 2 F\{a, b\a+ 1; gz)] 


(i a)^ _ ^ z b p or || g 341 1 it i s a i so US ed 


the transformation formula 2 A (a, a; a 4 - 1; -~3j) = (1 — z) a 2 A (a. l;a + 1; 2 ) [38 


= az 
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condition has been implemented by adopting the Weitzenbock connection whose coordinate components are a function 
of the vierbein, so our torsion is the torsion of the Weitzenbock connection and the dynamical object of the theory 
is solely the vierbein. The theory consists of the Einstein gravity in its teleparallel representation supplemented by 
a minimally coupled scalar field with potential. The novel additional term of the theory is a non-minimal derivative 
coupling of the scalar field with the torsion scalar (which is a particular quadratic combination of the torsion that 
provides under variation with respect to the vierbein the Einstein tensor). The equations of motion are second order 
differential equations and the theory (action and field equations) is both diffeomorphism and local Lorentz invariant, 
while after adopting the Weitzenbock connection in order to proceed, the local Lorentz invariance is abolished. 

The master equation found in [27], determining the four-dimensional spherically symmetric solutions, has here been 
elaborated and appropriately approximated in order to find large-distance solutions. The solutions found here are 
non-linearized, so they are valid at any distance larger than some length scale defined by the parameters/integration 
constants of the problem (the region of the parameters/integration constants is not fine-tuned or particularly narrow). 
A dynamical systems analysis has been performed and elucidates in the space of dynamical variables describing all 
the local solutions the regions where the large-distance solutions inhabit. This offers an intuition for the generality of 
the approximation method applied in order to obtain the general solutions. 

Special solutions have been found, among which black holes also, when the integration constant of the scalar field 
takes a particular value. Asymptotically, these solutions get an AdS form. The corresponding scalar field is a secondary 
hair. There are branches where the asymptotic value of the potential coincides with the effective cosmological constant 
of AdS; also there are other branches where the cosmological constant coming from the potential is positive (therefore 
it is different than the effective cosmological constant of AdS) and this is due to the growth of a non trivial profile of 
the scalar field at infinity after interacting with torsion. 

Probably the most interesting solutions found are branches of general (from the point of view of the number of 
integration constants) spherically symmetric solutions. These solutions are attracted by stable fixed points at the 
asymptotic region of the phase portrait of the dynamical system. Among them, there is the general black hole solution 
described by the metrics (16.271) . (16.301) with their asymptotic forms (16.281) . (16.321) . The form (16.301) clearly shows that 
the solution is a hyperscaling violating black hole with positive hyperscaling violation exponent and dynamical critical 
exponent. Note that the solution obtained here is a general solution of our theory, while the majority of the Lifshitz 
or hyperscaling violating solutions of differing theories in the literature are special. Note also that the topology of 
the horizon here is spherical, while most of the existing Lifshitz or hyperscaling violating black holes have planar 
horizon topology. Due to the spherical symmetry the solution found may have astrophysical significance, e.g. may 
possess extra attraction at large distances. Actually, we have found for a freely falling photon in the static isotropic 
gravitational field an increase of its deflection compared to the Newtonian deflection. The scalar field <j) accompanying 
the solution is a primary hair which diverges at the horizon and is finite at infinity. Notice, however, that in all the 
black hole solutions found in this work, although the scalar field diverges at the location of the horizon, its energy 
density and the pressures are finite there. 
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